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ABSTRACT 

Continuous-time, random-walk Monte Carlo simulations of H2 formation on grains 
have been performed for surfaces that are stochastically heated by photons. We have 
assumed diffuse cloud conditions and used a variety of grains of varying roughness and 
size based on olivine. The simulations were performed at different optical depths. We 
confirmed that small grains (r ^ 0.02 /im) have low modal temperatures with strong 
fluctuations, which have a large effect on the efficiency of the formation of molecular 
hydrogen. The grain size distribution highly favours small grains and therefore H2 
formation on these particles makes a large contribution to the overall formation rate 
for all but the roughest surfaces. We find that at Ay — only the roughest surfaces 
can produce the required amount of molecular hydrogen, but by Ay — 1, smoother 
surfaces are possible alternatives. Use of a larger value for the evaporation energy of 
atomic hydrogen, but one still consistent with experiment, allows smoother surfaces 
to produce more H2. 

Key words: ISM: molecules - molecular processes 



INTRODUCTION 



Although molecular hydrogen is abundant in diffuse and 
dense regions of the neutral interstellar medium (ISM), the 
detailed mechanism of its formation is still not completely 
clear. At the low gas-phase temperatures typical of the inter- 
stellar medium, a gas-phase formation is not possible, and 
H2 is therefore formed on the sur faces of grains. An inte r- 
pretation of experime ntal studies jPirronello et alJll997al lbl 
iPirronello et alj 119991) indicates that the surface tempera- 
ture range over which efficient H2 formation occurs is very 
sma ll for olivine (6-10 K) and for amorphous carbon (13-17 
K) iKatz et al.lll999l) under diffuse-cloud conditions. The 
grain temperature for midsize grains in unshielded regions 
is however around 20 K. This means that molecular hydro- 
gen cannot be formed efficiently on such grains assuming 
that the olivine surfaces used in the laboratory are a re- 
alistic repr e senta tion of interstellar grains. The analysis of 
iKatz et alJ il999l) considered, however, only a single barrier 
for diffusion between binding sites and a single low binding 
energy betwee n physisorbed H atoms and the surface. In pre - 
vious papers llChang et alJl2005t ICuppen fc Herbsd l2005l) . 
we showed that the introduction of sites of different binding 
energy leads to an increase of the temperature range for ef- 
ficient H2 formation. The different sites are the result of a 
difference in local environment either due to the amorphous 



character of the grain or to the topological structure of the 
grain; i.e., surface ro ughness. These resu l ts are in agreement 
with the findings of ICazaux &: TielensI (j200^, who intro- 
duced extra chemisorption sites and made a distinction be- 
tween H and D atoms. They found a larger temperature 
range depending on the width and the height of the barrier 
between sites of physisorption and chemisorption. 



In this paper, we will study the formation of molec- 
ular hydrogen using similar continuous-time, random- 
walk Monte Carlo simulations as in our pr evious papers 
JChang et al.l l2005l : ICuppen fc Herbsd l2005l) while grains 
are heated by the interstellar radiation field. The photons 
hitting a grain give it a short heat impulse, resulting in 
temperature fluctuations especially for small grains. These 
temperature fluctuations have receiv ed considerable atten- 
tion in the literature (GrccnbcrE & Ho na urce lil976l : 
lAannestad fc Kenvon, 1979; Tabak. 19871: loraine fc Ll200lT) . 
where they are discussed with a varying degree of detail. 
While most of the studies simply use the heat capacity 
to calculate the temperature changes, Draine & Li (200j|) 
looked at the heating and cooling of small silicon clusters 
and PAHs by considering the infrared emission of the grains, 
taking into account the infrared bands of the grain material. 
Since our main goal is to study molecular hydrogen forma- 
tion, we will use a relatively simple model for the heating 
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and cooling. As we will show, this model is capable of re- 
producing the detailed method to a large extent. 

Since both the heating and the recombination of H 
atoms into H2 are stochastic processes, the Monte Carlo 
technique is an ideal method to study how both processes 
influence each other. During the simulations, the individ- 
ual H atoms are followed as they land on the grain surface, 
undergo their random walk, and evaporate. If two atoms ar- 
rive at the same position on the grain, they react to form 
molecular hydrogen. As these processes occur, photons are 
hitting the grain, leading to a pulsed temperature increase 
and subsequent decrease for the smaller grains. The radia- 
tion field and the absorption and emission coefficients used 
in the simulations are explained in Section |5| The times at 
which the possible events - hopping, evaporation, and de- 
position of If atoms, and heating of the grain - occur are 
randomly determined using the rates of the processes and a 
random number . The details of the Monte Carlo procedure 
are discussed in IChang et all ll2005l) . We briefly summarise 
them in Section |21 which also explains the heating and cool- 
ing of the grain in more detail. Section|l]contains our results 
for the fluctuations in temperature, while Section |21 shows 
their effect on the formation efficiency of II2 as well as the 
effects of grain temperature vs. grain size and immediate 
grain desorption following formation of H2. Our results are 
put into the context of the ISM in Section |S| Finally, our 
conclusions are discussed in Section |7| 



2 THE PHOTON FLUX 

The rate of photons absorbed by a grain depends on its po- 
sition in the cloud and on its optical properties. We assume 
the distribution of photons absorbed by the grain per sec- 
ond per wavelength interval as a function of wavelength to 
be given by 



Pa = vrr IxQabs{\)Dx, 



(1) 



with r the radius of the grain, Ix the interstellar radiation 
field in photons per unit area per second per unit wave- 
length, Qabs the wavelength-dependent absorption coeffi- 
cient, and Dx the reduction factor due to attenuation of 
the radiation field by dust. 

For the interstellar radiation field Ix we use different 
expressions in different regimes. In the high photon energy 
regime, wi t h wav elengths 91.2-250 nm, the approximation of 
ISternber3 il988t) is used. The lower photon energy regime, 
with wavelengths 250 nm-1 cm, is divided into f our different 
regions using the expressions from lZucconi et'a l. (2001). 

At low photon energies, the absorption coefficient, Qabs, 
is approximated using 



Qabs(-\) ~ — — Im 
A 



(2) 



jBohren fc HuffmanlHoil . where m is the complex refrac- 
tive index: 



m = n + ki. 



(3) 



For the optical cons tants n and fc, we used tabulated val- 
ues for olivine from lOssenkopf et al.1 (|l992l ). Since Eq. ^ 
is only valid for 2ll <g 1 we only used this approximation 
for Qabs for A > 1000 nm. For smaller wavelengths, the r 



dependence of Qabs becomes different. Here we used a set 
of empirical expressions to approximate the absorption co- 
efficients found bv .Draine fc Lo c (1984) and Dr aine 
They showed the absorption coefficients for silicates and car- 
bonaceous material to be significantly different over a large 
energy range. We only used the silicate results here. 

For t he shielding of the radiation field, we followed 
iDraine ~Bertoldi, (.199(i) . who derived the attenuation fac- 
tor, Dx: 



Dx = exp -3.7 



Ax 



A, 



(4) 



based on the attenuation at 1000 A. Since we preferred to use 
Ay in stead of A^oooa we took ^1°°°^ = 4.70 from I Whitteil 
i2003h to arrive at the alternative expression 

Ax 



Dx = exp f-0.8-^A^ 
V Av 



(5) 



We then used the table of -j^ values in iMathij (ll99Cl) and 
IWhitt et (2003). 

The left panel of FigureQgives the total fiux of absorbed 
photons Px in photons per second for grains of two different 
sizes and at different visual extinctions as a function of en- 
ergy. The shape of the curves for different grain sizes is the 
same for low energies because of the constant dependence 
stemming from the linear dependence of Qabs (A). This re- 
lationship no longer holds in the high energy regime where 
the curve shapes are different for different grain sizes. This 
regime can be divided in two ranges: Px depends on more 
than from approximately 1-7 eV while the dependency is 
less strong for energies higher than 7 eV. The graph further 
shows that dust attenuation mainly affects the high energy 
regime. 



3 THE SURFACES AND THE MODEL 



In our previous paper (ICuppen fc Herbstl '2005'). we per- 
formed our Monte Carlo simulations of the recombination of 
hydrogen on different surfaces of varying surface roughness. 
These surfaces were generated using a separate Monte Carlo 
simulation program starting from a ffat surface with a square 
grid to represent binding sites on a grain. The resulting sur- 
faces were characterised by a matrix indicating the topo- 
logical roughness of each lattice grid. We used four different 
surfaces: surface (a) is completely ffat while surfaces (b), (c), 
and (d) have protrusions and "holes". In addition, surface 
(c) has several features known as islands and surface (d) 
is very rough with height differences of several monolayers 
and no distinct shape for the protrusions. If a hydrogen atom 
lands on the irregular surface, its local hopping and evapora- 
tion energy are calculated based on its number of horizontal 
'grain' neighbours, i (see Fig. 2 in lCuppen fc HerbstI (j 2005D . 
For every horizontal neighbour, there is an extra interaction 
energy El- This method effectively results in surfaces with 
a maximum of five different binding sites. The present paper 
uses the same concept to construct the three surfaces used 
in the new simulations. Surface I has only one type of site 
and is therefore ffat. Surface II has a roug hness comparable 
to surface (c) in lCuppen fc Herbsd (|200^, but we only dis- 
tinguish two types of sites. The first type of site has zero 
(i = 0) horizontal neighbours while the second type, which 
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Figure 1. (left) The absorbed photon flux in photons per second per wavelength interval and (right) the relation between the random 
number and the photon energy for different grain sizes, r = 0.01 and 0.1 /xm, and at different visual extinctions, Ay = 0, 0.5 and 1 (see 
legend, which applies to both panels). 



can have one or more horizontal neighbours according to the 
method used to generate the surface, is still assumed to have 
only o ne neighbour (i = 1 ). Sur face III derives from surface 
(d) in lCuppen fc Herbsd ll2005l) and has five different types 
of sites, with to 4 horizontal neighbours. 

Figure |21 shows representations of the non-flat sur- 
faces II and III. Note that the colours used in this flgure 
indicate the different si tes, in contrast with Figure 1 in 
ICuppen fc Herbs^ ^2005^ where the colour coding is used 
to indicate the different heights. For surface II, the black 
sites correspond to i = and the white sites to i = 1, whilst 
for Surface III, the lighter the colour, the greater the number 
of horizontal neighbours i. The hopping (_Et.) and evapora- 
tion {Ed) energies of these sites, labelled by the number of 
neighbours i, are given by 

Eb,r = Et + iEL, (6) 

ED,r = Ed+iEl, (7) 

respectively. The energies Ei,, Ed, and El used here for H 
and H2 are given in Tabled for the different surfaces along 
with the total number of different sites, itot and the dis- 
tribution of these different sites n^. The lateral bond. El, 
is chosen to be 40% of the evaporation barrier, Ed- Our 
previous paper dCuppen fc Hcrbst 2005) showed that the 
recombination efficiency for surface III is close to unity for 
temperatures around 20 K. With this choice of the lateral 
bond, we can distinguish three cases: zero, low, and high effi- 
ciency for the surface temperatures of interest in the absence 
of fluctuations. 

Each Monte Carlo cycle starts by checking which event 
will occur next. This is done by comparing the times at 
which all future events occur; heating of the grain, cool- 
ing of the grain, deposition of an H atom, or hopping or 
evaporation of a particle on the surface. All these events 
and their corresponding times are discussed in the following 
three sections. We start by explaining the events involved in 
the recombination of molecular hydrogen while Sections 13. 21 
and 13.31 contain discussions of the heating and the cooling 
of the grain, respectively. 




Figure 2. The distribution of the olivine surface sites for (top) 
itot = 2 and (bottom) itot = 5. The black area indicates the sites 
with i = and the lighter colours are for increasing i. 



3.1 Desorption, hopping, and evaporation of H 
and H2 

The deposition rate of H atoms (s~^ site~^) is calculated 
from the gas abundance using 
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Table 1. The parameters characterising the different surfaces 
used. 
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0.126 
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0.024 



1 taken from lKatz et all 111 

0.77gn.Ho. where .77 is the ratio between h 



3l200(t . 



and Er 



: 0.4Er 



"^ iBiham et aO J2i 

wn(H) 



Rdep — 



4p 



(8) 



where n(H) is the concentration of atomic hydrogen, p is the 
surface site density, and v is the average velocity of the H 
atoms in the gas: 



(9) 



with Tgas the temperature of the gas and m the atomic mass. 
We us e n(H) = 102 cm' ^ Tgas = 60 K, and p = 2 x lO" 
cm "2 iBiham et alJl200j) . The time between two deposition 
events is determined using a random number X between 
and 1, and the deposition rate: 



dep 



RdepNs 



(10) 



where A^s is the number of sites on the grain surface. 

When a deposition event occurs, the deposition site is 
determined using another random number. If the atom lands 
on a bare site of the grain, it simply lands there and the time 
at which the particle will do its next event is determined. 
This is done in a way similar to determining Atdep'- 



ievent 



ln(X) 



(11) 



Rhop ^" Reva 

where t is the current time in the simulation and Rhop and 
Reva are respectively, the hopping rate and the evaporation 
rate of the particle. These two rates are given by 



R 



hop,i = f CXp 



and 



Reva,i = Z^exp 



Eb,i 

Ed. 



{'¥) 



(12) 



(13) 



respectively, where T is the temperature of the dust and u 
is a trial frequency here chosen to be lOi^ s~i (Table0. A 
second random number is used to determine which event, 
hopping or evaporation, will occur at that time. 



If the particle (H) lands on top of another particle (H or 
H2), it is checked if a reaction between the landing species 
and species already there can occur. If so, the old particle is 
replaced by the new molecule, otherwise the incoming parti- 
cle remains on top of the old one with hopping and evapora- 
tion rates, Rhop,o and Reva,o, respectively. In both cases, the 
time of the next event for the new or newly formed parti- 
cle is determined. We assume initially that all H2 molecules 
that are formed stay on the surface upon reaction and later 
evaporate. T his corresponds to the parameter choice n — 1 
introduced bv lKatz et alJlll999t) . where p, indicates the frac- 
tion of molecules that remains on the surface immediately 
following reaction. 

3.2 Heating of the grain 

We consider photons in the wavelength range 91.2 nm-lO 
mm. The rate (s~i) of photons absorbed by a grain is 



R 



'photon — 



PxdX 



(14) 



from which the time between two photon hits can be deter- 
mined in a manner similar to Eqs. IIOII and lllll . If the grain 
absorbs a photon, the energy E of the photon is picked from 
the distribution in Eq. Q using a random number. The right 
panel in Figure Q indicates the relation between the random 
number, X\, and the photon energy. This figure is obtained 
using the relation 



X, 



J9: 



PxdX 



(15) 



Rphoton 

which gives a value between and 1 for a given wavelength 
A. This value is difi'erent for different grain sizes and dif- 
ferent visual extinctions. The r dependence is due to the r 
dependence in the absorption coefficient. 

The new temperature is then obtained by solving 

E= / c{T)AT (16) 



for Tnew , 

to [ Aa nncstad 
iPurce 



where c(T) is the heat capacity . Acc ording 
iDraine fc Lil \ 



KenvoiJ 



and 



PurcelL 1,1976.') , olivine can be approximated by a Debye solid 
with a Debye temperature of 500 K. The low temperature 
heat capacity in the Debye model is 

c{T^O) = ^Nk(^Y , (17) 



with A*' the number of atoms in the grain. If we assume 
a density of 3.32 g cm"'' and a molecular mass of 153.3 g 
mol~^, the heat capacity in eV K~i becomes 

c(r) = 61.38r^^^ (18) 

with r in /im. 

Every time the grain temperature changes, the time of 
the next event for each particle on the surface needs to be 
adjusted from to Cuent according to the new temper- 

ature. This procedure is done using the equation 

RhopiTold) + RevaiTold) 



,new /lOld ,\ 

''event \''event 



Rhop(.Tnew^ ~\~ Revai^new^ 

Since the ratio between Rhop and Reva also changes, the type 
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Figure 3. The temperature fluctuations as a function of time for 
different grain sizes and Ay = 0. 



3.3 Cooling of the grain 

The temperature of the grain is recalculated at certain time 
intervals using 



At = 



c{T) 



dT, 



where ^ is given by 



"37 

—— — 4tv { IT a 
dt ^ ' 



(A) Bx (T) dA 



(20) 



(21) 



with Qem the emission coefficient and B the Planck function. 
We found that A.t{s) — max (lOr^ //Jphoton, 30) with r in 
/xm, is a good time interval, because it takes into account the 
difference in photon flux and the temperature fluctuations. 

We integrate Eq. II2H numerically over the range 91.2 
nm < A ^ 10 mm using Qem. (A) = Qabs (A) and solve 
Eq. 1201 to obtain the new temperature. 



4 TEMPERATURE FLUCTUATIONS 

Before studying the formation of molecular hydrogen, we 
first discuss the influence of the incoming photons on the 
grain temperature. Figure |3 shows how the grain tempera- 
ture fluctuates as a function of time for different grain sizes 
at Av = 0. It clearly shows that although smaller grains 
have fewer incoming photons per given time interval, the 
temperature fluctuations are much larger, which is a con- 
sequence of the term in the heat capacity. For a grain 
size (radius) of 0.005 /xm, fluctuations of up to 30 K occur, 
while for a grain size of 0.05 /xm, an almost constant temper- 
ature is obtained. Furthermore, it can be noticed that the 
minimum temperature reached before each photon hit is al- 
most independent of the height of the temperature increase 
given the rapid initial cooling. [Praine ( 200_3,) shows similar 
graphs for carbonaceous grains, in which the temperature for 
the grains has similar fluctuations both in amplitude and in 
modal temperature. 

To obtain the modal grain temperature and the ampli- 
tude of the fluctuations, we determined the temperature at 




20 

Temperature (K) 



Figure 4. The temperature distribution over time for Ay = 
and r = 0.01 ^m. The lines indicate the modal temperature 
(solid) and the two 99% levels (dashed). See the text for more 
detail. 




0.01 0.1 
r(nm) 

Figure 5. Modal temperature of the grain as a function of grain 
size for different visual extinctions. The bars indicate 99% confi- 
dence levels. 



flxed intervals chosen such that there were several photon 
hits. Figure 0] shows a histogram of the temperature distri- 
bution so obtained for Ay = and r = 0.01 /xm. The his- 
togram clearly has a maximum, at which the temperature of 
the grain occurs most frequently. This mode is indicated by 
a solid line, which is left of the centre of the distribution as 
could be expected considering the shape of the temperature 
vs. time peaks. Figure |^ gives these modal grain temper- 
atures as a function of grain size for Ay — 0, 0.5 and 1. 
The bars indicate the 99% boundaries; i. e. the tempera- 
tures between which 99% of the determinations fall. These 
boundaries are indicated in Figure |1] by dashed lines and 
are a measure of the amplitude of the fluctuations. Figure 
|K|clearly shows that the temperature fluctuations are larger 
for smaller grain sizes. All three temperature proflles show 
a maximum, which is at rm ~ 0.02 fj,Tn for Ay — and 
gradually moves to larger r for higher Ay. The maximum 
is sharper at Ay ^ 0, which indicates that the radiation 
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Figure 6. Modal temperature of the grain as a function of Ay 
for different grain sizes. 



causing this peak is in the UV since this region is strongly 
shielded at higher extinction. The effect is due to the r de- 
pendence of the absorption coefficient in the high energy 
ran ge, as discussed in S ection |21 

iLi fc Draind ||200 show a similar graph, where the av- 
erage temperatures of both graphite and silicate grains are 
plotted for different radiation field intensities. Comparison 
between the xmmp = 1 curve from that graph and the upper 
left panel of Figure 13 shows that the shape of both curves is 
the same, but that the temperatures in Figure|S]are slightly 
higher (1 K). This is a good agreement considering the differ- 
ence in detail used in the methods. Note that in the limited 
regime depicted bviTi fc Draine (2o63) (r > 0.01 fim) the 
modal and average temperatures coincide. 

Figure |H| gives the modal grain temperature but now as 
a function of Ay for three different grain sizes. The three 
curves show very different behaviour. The modal temper- 
ature of the smallest grain remains constant as a function 
of extinction but the amplitude of the fluctuations clearly 
decreases. The intermediate grain of radius 0.02 fim has a 
gradually decreasing modal temperature but the amplitude 
remains constant. It is interesting how the position of the 
mode in the distribution changes from the centre for Av = 
to the lower temperatures for Av — 1.5. Finally, the 0.1 
^m grain has almost no temperature fiuctuations and the 
curve shows a gradual decrease of the temperature for grains 
deeper into the cloud. Notice the crossing of the curves for 
0.02 fim and 0.1 fim around Av ~ 1, which is a consequence 
of the gradual moving of to larger values for higher visual 
extinction. 



5 EFFECT OF THE TEMPERATURE ON THE 
EFFICIENCY 

The previous section showed that not only the amplitude 
of the temperature fluctuations is grain-size dependent but 
also the modal temperature. This means that we should not 
limit our study to sizes where fluctuations play a role but 
study a wider range of grain sizes in order to obtain the 
complete size dependence of the recombination efficiency. 



o 0.4 



0.2 




10 12 14 16 
Temperature (K) 



20 22 24 



Figure 7. The recombination efficiency as a function of tempera- 
ture for the three different surfaces in the absence of fiuctuations. 
The closed markers all indicate simulations with /i = 1, while 
for Surface III, additional simulations are run and designated by 
open markers for fi = (upper curve) and fj, = 0.3 (lower curve). 



The definition of the recombination efficiency, rj, 
in this paper is 

2iVH, 



V 



used 



(22) 



where A''h is the number of hydrogen atoms that approach 
the surface in a given time interval and A'h2 is the number 
of molecules that leave the surface in that interval. Usually, 
rj is defined as 

2i?H, 



V 



(23) 



where is the rate with which H2 comes off the surface 
and F is the flux of incoming hydrogen atoms. Since the 
Monte Carlo method uses discrete numbers and not rates, 
as in rate equation methods, the evaporation rate cannot be 
determined directly. We therefore use Eq. 1221 which gives 
the same results as Eq. 12311 in steady-state conditions. 

We first consider the efficiencies as a function of the 
temperature for the three different surfaces introduced in 
Section |3| The results are shown in Figure Q and were ob- 
tained with a simulation in which temperature fluctuations 
were omitted. The size of the grains is large enough that the 
efficiency is roughly independent of it ijChang et alJl20oil) . a 
regime known as the accretion limit. These graphs will help 
us to interpret the results we obtain from the simulations 
with temperature fluctuations. The efficiency curve for sur- 
face I, flat olivine, has one narrow peak at 8 K, surface II 
(olivine with islands) has one broad peak (7-11 K) punctu- 
ated by one dip and surface III (very rough surface) has a 
very broad peak (6-18 K) with three dips. These dips are 
due to hydrogen molecules, which occupy the high energy 
sites (i > 0) and prevent hydrogen atoms from using these 
sites to make new molecules. At slightly higher temperatures 
the formed molecules evaporate and the sites are available 
again. Since we use = 1 all molecules remain on the sur- 
face aft er formation until they evaporate. In our previous 
papers dChang et 8111200,4 IciiDPen fc Herbstll200,'^^ we used 
H = and only one peak was found in all cases. The open 
markers in Figure |7| give simulation results for /i = and 
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Figure 8. The recombination efficiency as a function of the grain 
size for different surfaces and visual extinctions. 



Figure 9. The recombination efficiency as a function of Ay for 
different surfaces and grain sizes. 



H = 0.3 on Surface III, where the n dependence is signifi- 
cant. The value of 0.3 is chosen since it is the value found by 
iKatz et alJ lll998t) for H2 formation on olivine. Notice that 
for certain temperatures the difference in efficiency can be 
as large as a factor of 10. These differences are however lim- 
ited to very small temperature ranges so that the value of /i 
only influences the average efficiency over the range within 
a much smaller factor between the two extremes fi = and 
H = 1. The efficiency curves for ^ = 0.3 and fJ. — 1 are much 
closer. We use n — 1 for now. 

Figure |S] shows the recombination efficiency as a func- 
tion of grain size with the temperature fluctuations depicted 
as in Figure |K| for the three different surfaces. With our as- 
sumed surface site density of 2 x 10^* cm~^, we use array 
sizes of 25 x 25, 50 x 50, and 100 x 100 to represent the 
surfaces of grains of radii r — 0.005, 0.01, and 0.02 fim, re- 
spectively. The efficiency at a given temperature for larger 
grains is in all cases independent of size. So we can use a 
100 X 100 array for such grains as long as the temperature 
fluctuations are headed correctly. 

Consider first the results for the smallest grains: r = 
0.005 jj,m. From Figures |3 and |S| we know that these grains 
have a modal grain temperature of 8 K. According to Figure 
[7|the efficiencies at this temperature are approximately 1.0, 
0.3 and 0.1 for surfaces I, II, and III, respectively, without 
temperature fluctuations. Figure |S] shows however that for 
these small grains the recombination efficiency is similar for 
all surfaces and rather low. The temperature fluctuations, 
which occur frequently compared with the H-atom deposi- 
tion interval, increase the evaporation rate which causes the 
average residence time of atoms on the surface to be reduced 
to such an extent that the efficiency is much lower and the 
atoms are no longer available when the grain cools down. As 
the grain size becomes somewhat larger, the modal temper- 
ature rises and the efficiency actually decreases for surfaces 
I and II, showing the rise in the modal temperature to out- 
weigh the decrease in fluctuations. 



If we look at the larger grains, r ^ 0.05 /im, the fluc- 
tuations become even less salient and the efficiencies are as 
could be expected combining Figures |S] and |7| A grain of 
size r ^ 0.1 /xm at Av = 1 for example has a temperature of 
14 K which would result in efficiencies of approximately 0, 
0, and 0.65 for surfaces I, II, and III, respectively, without 
temperature fluctuations. Figure |S| shows that this is indeed 
the case. 

Figure |5] shows the efficiency again, but now as a func- 
tion of visual extinction for three grain sizes, and should 
be studied in combination with Figure |S] All graphs show 
an increase for increasing Ay, except the 0.1 /im curve for 
surface III. These increases in efficiency are due to the tem- 
perature decrease as a function of visual extinction for grains 
of 0.02 /xm and due to a decrease in fluctuations for grains of 
0.005 fim. For grains of 0.02 fim with surface III, the fluctu- 
ations gradually move inside the temperature range where 
a non-zero efficiency is obtained. This occurs at Ay ~ 0.7, 
which is the reason that for higher extinction the efficiency 
hardly changes. The minimum in the 0.1 pm curve for sur- 
face III corresponds to the dip in the efficiency curve (Figure 
|7|l around 15 K, since the temperature decreases from 17.5 
K to 13 K going from Av = to 1.5. 

Figure [TI71 shows the efficiency as a function of Ay for 
^ = and 0.3 for surface III; it can be compared with the 
lower left panel of Figure]^ which gives results for jj, = 1. 
We have showed already that the value of fi can have a 
strong effect on the efficiency at specific small temperature 
ranges, but have also mentioned that on average the effect 
is much smaller. The low efficiency ranges with /i = 1 are 
around 8 and 11 K, at which newly formed H2 molecules 
interfere with subsequent reactions. Although the grains of 
0.005 /xm have a modal temperature of 8 K, the fluctuations 
cause H2 to evaporate easily, resulting in a low effect of /x on 
the efficiency. For the other two grain sizes we see a visible 
but still small effect due to variations in /i. This could be 
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Figure 10. The efBciency as a function of the visual extinction 
for different values of fj, for surface III and different grain sizes. 

expected considering that the modal temperatures (12-20 
K) do not fall in the low efBciency dips. 



6 CONSEQUENCES FOR THE DIFFUSE ISM 

In the previous sections, we showed that the amplitude of 
the temperature fluctuations and the modal temperature of 
grains are strongly size dependent, and that both quantities 
have a strong effect on the efficiency of molecular hydrogen 
formation on olivine. While for the large grains the type of 
surface also has a strong effect on the efficiency, the rough- 
ness of the surface seems to have little effect for the smallest 
grains, since for those grains there is enough thermal energy 
for atoms to hop and evaporate even from the high energy 
sites because of the strong temperature fluctuations. 

This section discusses the consequences of these findings 
for molecular hydrogen formation in the diffuse interstellar 
medium. Usually, formation rates for molecular hydrogen 
are based on the assumption that all grains have a standard 
size of 0.1 ^m. In reality, there is a wide range of differ- 
ent sizes in the ISM. The distribution of these sizes highly 
favo urs small grains an d can be approximated by a power 
law llMathis et alJll977ll : 

ng(r)=cr-^-'^, (24) 

where Ug (r)dr is the conc entration of grains with radii in the 
range (r, r+dr in cm'^'t. lbraine fc Lee (1984) found a value 
c — 1.51 X 10 '^uh cm 2 for the proportion ality constant 
using rmin ~ 0.005 fim and rmax = 0.25 (iMa this et alJ 
|l2Z3) for the minimum and maximum radii. The gas-to-dust 
number ratio is found by integrating Eq. I24II over the range 
and dividing by uh to be 2.91 x 10' , which is significantly 
lower than the usually-assumed value of 10^ , a number that 
only considers 0.1 fim grains and therefore neglects the large 
portion of small grains. A somewhat hi gher value of 6.0 x 10® 
was obtained bv lLipshtat et aP (|20o3) who only considered 
grains down to 0.04 /xm in size. 

From the grain size distribution we can calculate the 
rate coefficient R in cm^ s~^ for molecular hydrogen forma- 
tion: 

X!""" ^ng{r)vri{r)^Tvr^dr 
R = , (25) 

where ^ is the sticking coefBcient of the atoms to the grain. 
For the low temperatures used here the sticking coefficient 



can be assumed to be close to one iBuch fc Czerminskll99ll : 
lAl-Halabi et alJl2002() . Eq. ^ can be simplified to 

where 

a = cJ / r){r)r ^'^dr. (27) 

Note that in the standard treatment where all grains are 
0.1 /^m in size and 77 = 1, q = 5 x 10~^^ cm^ s~^. This is 
in agreement with the estimate oi R — 10~ ^^ cm^ s~^ for 
diffuse-cloud conditions given bv I Jural lll974ll . which corre- 
sponds with a fa 2 X lO"^'^ cm^ s~^. 

Figure 1111 shows T]{r)r~^'^ as a function of the grain 
size using the data from Figure |S] For surfaces I and II 
the maxima for the efficiency and r~^'^ coincide at rmin, 
which obviously leads to a maximum for the smallest grain 
sizes used. Surface III, however, has its largest efficiencies for 
the largest grains. This results in a minimum for 'i]{r)r~^'^ 
around r = 0.01 fim and a maximum at intermediate grain 
sizes, depending on Av- Figure [TTI clearlv shows that the 
higher efficiency for large grains is of almost no importance. 
For flat grains and those with a minimal roughness, the 
overall H2 formation is mainly determined on grains with 
r ^ 0.01 fim, while for grains with a high roughness, the 
intermediate range of 0.01 ^ r ^ 0.08 /im is most impor- 
tant. Since surfaces I and II have their maximum value for 
r;(r)r~^'^ at rmin, the choice for the lower limit in grain size 
distribution influences the value for a. Generally the min- 
imu m size for olivine grains is considered to be r = 0.005 
Hm jMathis et"al]ll977l: iDraine fc Le^ll984^ ■ which is the 
value we take. lL^^Dr£dn?l2oSh concluded, however, that 
up to ~ 10% of interstellar Si could be in r ^ 15 A sili- 
cate grains. We did not consider these ultra-small particles 
in our calculations because, in our view, such particles need 
to be treated as a special case for several reasons: (i) they 
consist of up to only a few hundred molecular units and 
therefore may not retain their bulk properties, such as the 
heat capacity that we use to calculate the heating and cool- 
ing; (ii) the small number of units on the surface means 
that our smooth and rough topologies are imprecise, (iii) 
the efficiency of st icking by weakly bound adsorbates is un- 
likely to be high jHerbst fc Dunbajll99l|) and it is likely 
that chemisorption is needed for high sticking efficiencies; 
and (iv) strong adsorbate-substrate bonds may be necessary 
for H atoms to withstand the very large temperature pulses 
that photons will cause. Such considerations will also be 
necessary to treat molecular hydrogen formation on PAH's. 

Table|21shows our numerically integrated a coefficients. 
Note that the efficiencies used for the integration are only 
determined for one gas temperature and one gas concentra- 
tion of H. We did not check for an additional dependence of 
the efficiency on these two quantities. Considering that the 
efficiencies used to obtain the values in Table |5| are much 
lower than unity, the a coefficients are surprisingly close to 
5 X 10"^^ cm'^ s~^. But we must remember that the large 
number of small grains gives a very large surface area for 
reaction. Indeed, if we consider a constant high efficiency 
of one over the whole grain size range, an a coefficient of 
1.53 X 10~^^ cm'^ s~^ is found. Of course our computed val- 
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Figure 11. The recombination efficiency X r~^'^ (or integrand 
in Eq. 1271 ') as a function of the grain size for different surfaces 
and visual extinctions. 



Table 2. a coefficient in cm^ s ^. 




Table 3. a coefficient in cm^ s ^ at a constant temperature. 
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Eo is increased by 35 K. 

ues are lower. From Table |5| we can see that only our rough 
surface III is able to produce a sufficient amount of molecular 
hydrogen at Av = 0, while the other surfaces are possible 
alternatives, with a a factor of six below Jura's value, at 
Av ~ 1, which pertains at more internal regions for diffuse 
clouds. This difference is caused by the intermediate grains, 
which have a very low efficiency for surface I and II, but 
have appreciable values for surface III. 

The low efficiency of surfaces I and II for intermediate 
grains is due ultimately to the low evaporation energy for H 
given bv lKatz et alj il999ft . In their paper they do indicate 
that their stated energy is more precisely only a lower limit 
within a range of 3 meV. Using a higher value results in 
larger a coefficients for surfaces I and II since H evaporates 
more slowly, and makes these surfaces more suitable candi- 
dates for efficient molecular hydrogen formation, especially 
at non-zero extinction. The last row in Table |5| gives the a 
coefficients for surface II if the value oi Ed is increased by 
the extra 3 meV (35 K, an increase of 9%.) We see that there 
is a clear increase in the a coefficient of a factor of a few for 
low Av that gets larger as the extinction increases, so that 
for Av > 0.5 the rate coefficient is high enough to produce 
a sufficient amount of molecular hydrogen in diffuse clouds. 

Table 121 gives values of the a coefficient for simulations 
performed at a constant temperature. We choose the modal 
temperatures of the standard grain of r = 0.1 ^m, which are 
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17.6, 15.3, and 13.9 K for Av = 0, 0.5 and 1, respectively. 
In this way, we can separate out the effects of the low modal 
temperature and stochastic heating on small grains. No H2 
was formed within our simulation time for the completely 
flat surface I and slightly rough surface II at Av = 0; we 
can therefore only give an upper limit for these cases. The 
higher values in the case of the stochastic heating (Table |2l 
stem from the lower modal temperature for the small grains. 
Surface III, on the other hand, has a small increase in the 
formation rate at constant temperature due to the higher 
efficiency for small grains in the absence of the fluctuation 
in temperature. Here the increase in modal temperature is 
not as important because of the large temperature range for 
high efficiency. 



7 CONCLUSIONS 

Previous calculations of the rate of molecular hydrogen for- 
mation on interstellar grain surfaces in diffuse clouds have 
ignored the subtle role played by radiation in determin- 
ing the surface temperature and its ffuctuations as func- 
tions of grain size . As shown by earlier investigators (e.g. 
iDraine fc Lil i200lll 'l and confirmed here, the temperature of 
grains fluctuates strongly for smaller grains (r < 0.02 ^m) 
but not for larger ones, and the modal temperature peaks 
for grains in the middle range of size (0.02 - 0.05 /xm). We 
have shown that the efficiency of formation of H2 depends 
both on the modal surface temperature and its fluctuations. 
According to the type of material and its smoothness or 
roughness, each surface has a range of temperature where 
the efficiency is high. For olivine grains, which are consid- 
ered here, a smooth surface has only a very small range of 
surface temperature for high efficiency 6-9 K - while a very 
rough surface, with half of its binding sites different from 
the norm, allows efficient H2 production up to about 20 
K. For small grains, photons pulse the temperature above 
the range of efficiency, thus reducing it, although the low 
modal temperature aids H2 production, especially for flat 
and moderately rough olivine grains. The net result for these 
two types of surfaces is an enhancement in the efficiency for 
smaller vs. larger grains although the overall efficiency re- 
mains low at all sizes. The relatively high efficiency for small 
grains is augmented by the high overall surface area given 
the distribution of grain sizes, so that the production of H2 
is dominated by the smaller grains. For very rough grains, 
on the other hand, the efficiency is greater for larger grains 
since fluctuations are minimal and the modal temperature is 
still within the range of high efficiency, while it is above this 
range for smoother surfaces. Here the countervailing larger 
surface area for smaller grains results in a peak H2 produc- 
tion for grains of intermediate sizes. 

A calculation of the rate coefficient R for hydrogen for- 
mation in diffuse clouds for which grains from a minimum 
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radius of 0.005 /im to a maximum radius of 0.25 /im are 
included show s that in ord er to equal or exceed the value of 
R deduced bv I.Txiral (|l974l ^ to reproduce the H2/H balance 
in diffuse clouds, one needs a very rough surface of olivine 
if Ay = 0. At a slightly higher visual extinction of 1.0, flat 
and moderately roug h sur f aces l ead to values of _R a factor 
of six lower than the I Jural (|l974') value. Since the latter vi- 
sual extinction is quite reasonable for internal portions of 
diffuse interstellar clouds, it appears that the rate of forma- 
tion of H2 on olivine grains of widely differing topologies can 
possibly account for the production of this species although 
the rougher surfaces are better candidates. This conclusion 
rests ultimately o n the energy need ed for evaporation of H 
atoms deduced bv lKatz et al.l lll999f) for an olivine surface; a 
somewhat higher value, also consistent with their data, im- 
proves the suitability of the smoother surfaces for efficient 
hydrogen production, especially if the extinction is greater 
than zero. If the effects of the radiation field are removed, 
on the other hand, the difference between the smoother and 
rougher surfaces remains large. 

In our calculations, we have only considered grains 
larger than 5 nm in radius consisting of olivine, and we 
have only treated a standard radiation field. The inclusion of 
smaller silicate particles is an obvious extension, as is the in- 
clusion of carbonaceous grains, which would also allow us to 
consider much smaller particles, including PAHs. For these 
very small particles, temperature fluctuations will be severe, 
and strong H atom-surface bonds, known as chemisorption 
JCazaux fc Ti elens 200^), must probably be considered to 
allow H atoms to remain on the grains at the peak temper- 
atures. It is possible, if barriers to chemisorption are low, 
that the inclusion of small silicate particles and PAHs will 
increase the calculated values of the rate coefficient R for 
diffuse clouds, and we intend to investigate this effect in the 
near future. On the other hand, an enhanced radiation field, 
as associated with photon-dominated regions, will increase 
the frequency of fiuctuations in temperature for small grains, 
and likely decrease the overall rate of H2 formation although 
the effect has yet to be studied. 
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